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ABSTRACT 

We present the Suzaku spectroscopic study of the Galactic middle-aged supernova remnant (SNR) 
IC 443. The X-ray spectrum in the 1.75-6.0 keV band is described by an optically-thin thermal 
plasma with the electron temperature of ~ 0.6 keV and several additional Lyman lines. We robustly 
detect, for the first time, strong radiative recombination continua (RRC) of H-like Si and S around at 
2.7 and 3.5 keV. The ionization temperatures of Si and S determined from the intensity ratios of the 
RRC to He-like Ka line are ~ 1.0 keV and ~ 1.2 keV, respectively. We thus find firm evidence for an 
extremely-overionized (recombining) plasma. As the origin of the overionization, a thermal conduction 
scenario argued in previous work is not favored in our new results. We propose that the highly-ionized 
gas were made at the initial phase of the SNR evolution in dense regions around a massive progenitor, 
and the low electron temperature is due to a rapid cooling by an adiabatic expansion. 
Subject headings: ISM: individual (IC 443) — supernova remnants — radiation mechanisms: thermal 
- X-rays: ISM 



1. INTRODUCTION 

IC 443 (G189. 1+3.0), a Galactic supernova remnant 
(SNR) at a distance of 1.5 kpc (Welsh & Sallmen 2003), 
is located near the Gem OBI association and a dense gi- 
ant molecular cloud (Cornett et al. 1977) with OH maser 
emission (Claussen et al. 1997). These facts strongly 
suggest that the remnant originated from a collapse of 
a massive progenitor. Braun & Strom (1986) proposed 
that the shock wave has expanded into the pre-existing 
wind-blown bubble shell. This was confirmed by a kine- 
matical study of the optical filaments (Meaburn et al. 
1990). A comprehensive X-ray study of IC 443 was first 
made with the Einstein and HE AO- A 2 satellites (Petre 
et al. 1988). They estimated the SNR age to be ~3000 yr. 
Recently, Troja et al. (2008) derived the age of -4000 yr 
from the morphologies of the shocked ejecta and inter- 
stellar medium (ISM) revealed by XMM-Newton. Thus, 
IC 443 is a middle-aged SNR. 

Using ASCA, Kawasaki et al. (2002) found that the 
ionization degrees of Si and S were significantly higher 
than those expected from the electron temperature of 
the bremsstrahlung continuum. Therefore, it was argued 
that the plasma is in overionization (kT z > kT e , where 
kT z and kT e are ionization and electron temperatures, 
respectively). On the other hand, Troja et al. (2008) 
found that the plasma is in collisional ionization equilib- 
rium (CIE), or the overionization is only marginal based 
on XMM-Newton data. These two controversial results 
hinge upon small differences, if any, between the esti- 
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mated electron and ionization temperatures. 

In this letter, we investigate whether the overionized 
plasma is really present or not, by utilizing the superior 
spectral capabilities for diffuse sources of X-ray Imag- 
ing Spectrometers (XIS: Koyama et al. 2007) aboard the 
Suzaku satellite (Mitsuda et al. 2007). If present, we will 
study the plasma condition quantitatively to discuss the 
possible origin of the overionization. 

2. OBSERVATION AND DATA REDUCTION 

The northern part of IC 443 was observed with Suzaku 
on 2007 March 6 (Observation ID = 501006010). Three 
XISs located at the foci of the independent X-Ray Tele- 
scopes (XRT: Serlemitsos et al. 2007) were operating 6 . 
Two of the XISs are Front-Illuminated (FI) CCDs and 
the other is a Back- Illuminated (BI) CCD. The XIS was 
operated in the normal full-frame clocking mode with a 
spaced-row charge injection technique (Uchiyama et al. 
2009) during the observation. For the data reduction 
and spectral analysis, we used the HEADAS software 
package of version 6.5 and XSPEC version 11.3.2, re- 
spectively. The XIS data of revision 2.0 were employed, 
but reprocessed using the xispi software and the latest 
calibration database released on 2009 February 3. Af- 
ter the screening with the standard criteria 7 , an effective 
exposure of —42 ksec was obtained. 

3. ANALYSIS AND RESULTS 

Figure Q] shows the vignetting-corrected XIS image in 
1.75-3.0 keV, the energy band including the major lines 
of K-shcll emissions from Si and S. We extracted the 
spectrum from the brightest region, the rectangular with 
an angular size of 7' x 5' (see Fig.[lJ. This approximately 
corresponds to the "Center" region of Kawasaki et al. 
(2002). The XIS spectra of the two FIs in the entire 
energy range were made by subtracting the non X-ray 

6 Suzaku carried four XISs, but one is now out of operation due 
to a damage possibly by an impact of a micro-meteorite 

7 http: / /heasarc. nasa.gov/docs / suzaku / proccssing/criteria_xis.html 
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Fig. 1. — Vignetting-corrected XIS image of the northern part 
of IC 443 in the 1.75-3.0 keV band, shown on an intensity color 
scale. The coordinates (R.A. and Dec.) refer to epoch J2000.0. 
The data from the three active XISs are combined and smoothed 
with a Gaussian kernel of a = 25". The yellow square and the 
white rectangle indicate the XIS field of view and the region used 
in our spectral analysis, respectively. The optical Digitized Sky 
Survey image is overplotted in contour. 
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Fig. 2. — Full-band XIS-FI spectrum where the NXB is sub- 
tracted. The energies of prominent Ka emission lines from specific 
elements are labeled in the panel. 



background (NXB) constructed with the xisnxbgen soft- 
ware. The spectra were merged to improve the statistics, 
because the response functions are almost identical each 
other. Figure [5] shows the resultant spectrum. We can 
see several prominent lines of Ka emission from He- and 
H-like ions (hereafter, Hea and Lya). The centroids of 
the Lya lines were measured with a Gaussian line model, 
and compared with the canonical values of the Astro- 
physical Plasma Emission Database (APED: Smith et al. 
2001). The averaged center energy difference was +4 eV. 
Therefore, we added a 4 eV offset in the FI spectrum. 
The BI spectrum was made with the same procedure as 
the FIs, but an offset of —10 eV was added to correct the 
energy scale. In order to examine the ionization states 
of Si and S, we hereafter focus on the spectrum in the 
energy range above 1.75 keV (Fig. [3]). Detailed studies 
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Fig. 3.— (a) XIS spectrum in the 1.75-6.0 keV band. Black 
and red represent FI and BI, respectively. Individual components 
of the best-fit model for the FI data are shown with solid colored 
lines: blue, green, and gray are the VAPEC, Gaussians (Si-Lya, Si- 
Ly/3, S-Lya, and Ar-Lyct), and CXB, respectively. The lower panel 
shows the residual from the best-fit model. Two hump-like features 
are clearly found around the energies of ~2.7 keV and ~3.5 keV. 
(b) Same spectrum as (a), but for a fit with RRC components of 
H-like Mg, Si, and S (magenta lines). The residuals seen in (a) are 
disappeared. 



including the lower energies will be reported in a separate 
paper. 

We first fitted the spectrum with a model of a thin- 
thermal plasma in CIE state (a VAPEC model). The 
abundances (Anders & Grevesse 1989) of Si, S, and Ar 
were free parameters, while the Ca abundance was tied 
to Ar. Interstellar extinction was fixed to a hydrogen 
column density of Ah = 7 x 10 21 cm~ 2 with the solar ele- 
mental abundances, following Kawasaki et al. (2002) and 
Troja et al. (2008). The cosmic X-ray background (CXB) 
spectrum was approximated by a power-law model with 
photon index of T = 1.412 and the surface brightness 
in the 2-10 keV band of 6.4 x 10~ 8 erg cm~ 2 s _1 sr _1 
(Kushino et al. 2002). Since IC 443 is located in the 
anti-Galactic center direction, contribution of the Galac- 
tic ridge X-ray emission was ignored. In the initial fit, we 
found a significant inconsistency between the FI and BI 
data around the energy of neutral Si K-edge (1.84 keV). 
This is due to the well-known calibration issue of the 
XIS. Since the calibration for the FI CCDs is currently 
far better than the BI, we decided to ignore the energy 
band below 1.9 keV in the BI spectrum. This fit leaved 
further large residuals, in both spectra, at the energies 
of S and Ar Lya lines, and hence was rejected with the 
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TABLE 1 

Best-fit spectral parameters 



Component 


Parameter 


Value 


CIE (VAPEC) 


kT e (keV) 




0.61 (0.59-0.64) 




Zg; (solar) 




0.82 (0.78-0.85) 




Zq (solar) 




1.7 (1.6-1.8) 




Z\ F (solar) 




2.5 (2.2-2.8) 




VEM a (10 12 


cm 5 ) 


6.4 (6.3-6.6) 


Additional components 






E (keV) b 


Flux c 


Line 


Si Lye? 


2.006 


2.8 (2.7-2.9) 




Si Ly/3 


2.377 


0.21 (0.14-0.29) 




S Lya 


2.623 


0.84 (0.79-0.90) 




Ar Lya 


3.323 


0.11 (0.085-0.14) 


RRC 


H-like Mg 


1.958 


1.2 (1.0-1.5) 




H-like Si 


2.666 


2.2 (2.0-2.3) 




H-like S 


3.482 


0.46 (0.41-0.51) 



Note. — The uncertainties in the parentheses are the 90% 
confidence range. "Volume emission measure, J n e n p dV j '(47r D 2 ) , 
where n e , n p , V, and D are the electron and proton densities 
(cm -3 ), the emitting volume (cm 3 ), and the distance to the source 
(cm), respectively. b The fixed energy values of the line cen- 
ter or the K-edge. c Total flux in the unit of 10 — 4 photon cm -2 s — 1 . 

X 2 /dof of 935/270. 

We, therefore, applied a model of one-temperature 
VAPEC (CIE plasma) plus narrow Gaussian lines at 
2006, 2623, and 3323 eV, the Lya energies of the Si, 
S, and Ar, hydrogenic ions, respectively. This process 
is essentially the same as Kawasaki et al. (2002). Ly/3 
lines of the same elements were also added, but they 
were not significant except for Si (2377 eV). The result 
is shown in Figured, with the best-fit CIE plasma tem- 
perature of kT e ~ 0.94 keV. Although the x 2 /dof value 
was significantly reduced to 657/266, the model was still 
unacceptable. In fact, apparent hump-like residuals are 
found around ~2.7 keV and ~3.5 keV, in addition to the 
systematic model excess in the continuum at the energies 
above ~4.5 keV. 

We checked systematic error due to the CXB fluctu- 
ation by allowing the CXB intensity as a free parame- 
ter. The fit, however, did not improve at all. We also 
tried two-component VAPEC models with independent 
temperatures and abundances, but these were rejected 
with x 2 /dof = 737/265. The results were essentially the 
same as figure the hump-like residuals remained at 
~2.7 keV and ~3.5 keV. No additional CIE component 
nor non-equilibrium ionization (NEI) plasma model re- 
moved the hump-like residuals with significant improve- 
ment of the x 2 value. 

At the energies of the humps, no emission line candi- 
date from an abundant element is found. However, the 
energies are consistent with the K-shell binding poten- 
tials (I z ) of the H-like Si (2666 eV) and S (3482 eV). 
Therefore, the humps are likely due to the free-bound 
transitions to the K-shell of the H-like Si and S. A for- 
mula for the spectrum of radiative recombination contin- 
uum (RRC) is found in equation (21) of Smith & Brick- 
house (2002) 8 . When the electron temperature is much 
lower than the K-edge energy (kT e <C I z ), this formula 
is approximated as; 

£(*.)««p(-^). *»*,>/.. a) 

8 http:/ /cxc. harvard.edu/atomdb/physics/plasma/plasma. html 
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Fig. 4. — Emissivity ratios of H-like RRC to He-like Ka lines 
as a function of ionization temperature (fcT z ), predicted by the 
plasma radiation code of Masai (1994) for recombining plasma with 
electron temperature (kT e ) of 0.6 keV. Black and red solid lines 
represent Si and S, respectively. For comparison, the same ratios 
for CIE (kT z = kT e ) plasma (the APEC model: Smith 2001) arc 
indicated by dashed lines. The horizontal dotted lines represent 
the 90% upper and lower limits of the observed values. 

Thus, the width of the RRC structure depends on the 
electron temperature. We added the RRC of Equation [1] 
for H-like 9 Si and S. The kT e values of the RRC were 
linked to that of the VAPEC component. Then, the fit 
was dramatically improved to an acceptable % 2 /dof of 
290/264. Although the hump- like residuals were com- 
pletely removed by this fit, we further added the RRC 
model of H-like Mg. This step is reasonable because Mg 
is more abundant than Si and S in the solar abundance 
ratio of Anders & Grevesse (1989) and the K-edge en- 
ergy of H-like Mg (I z — 1958 eV) falls into the analyzed 
band. The x 2 /dof value was significantly reduced to 
267/263, which gives an F-test probability of ~ 3 x 10~ 6 . 
The best- fit parameters and model are given in Table [1] 
and Figure [3Jd, respectively. The spectrum was also fit- 
ted with the independent electron temperatures for the 
VAPEC and RRC components, but the values are consis- 
tent with each other within their statistical uncertainties. 

4. DISCUSSION AND CONCLUSION 

We have found that the 1.75-6.0 keV spectrum cannot 
be represented with CIE nor NEI plasma alone, but need 
the additional fluxes of the Lyman lines and the H-like 
RRC of Si and S (and possibly Mg). This is the first 
detection of clear RRC emissions in an SNR. In the fol- 
lowing, we quantitatively discuss the implications of our 
spectral results. 

From the best-fit model in Table [1] the flux ratios of 
H-like RRC to Hc-like Ka (Hea) line (Frrc I Ffiea) are 
given to be 0.28 (0.24-0.30) and 0.18 (0.15-0.20) for Si 
and S, respectively. These are compared, in Figure 21 
with the modeled emissivity ratios by the plasma ra- 
diation code of Masai (1994) for the electron tempera- 
ture of 0.6 keV. We find that the large observed ratios of 
-Frrc / -FHea are significantly above those in the CIE case 
(kT z = kT e ), but can be reproduced in the overioniza- 
tion case (kT z > kT e ). The ionization temperatures of 
Si and S are determined to be ~ 1.0 keV and ~ 1.2 keV, 

9 "H-like RRC" refers to the free-bound emission due to electron 
captures by fully-ionized ions into ground state of H-like ions 
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respectively. This is, therefore, the firm evidence of the 
overionized (recombining) plasma. 

The overionization claim for IC 443 was first argued 
by Kawasaki et al. (2002). Using ASCA data, they de- 
rived kT z to be ~ 1.5 keV from the S Lya/Hea flux 
ratio compared with the predicted emissivity ratio in the 
CIE plasma code. Troja et al. (2008) adopted the same 
analysis procedure to the XMM-Newton spectrum, and 
claimed that kT z obtained from the line flux ratio was 
nearly same as the bremsstrahlung temperature (kT e ). 
However, since the RRC process is accompanied with 
electron captures to the excited levels, as given in Equa- 
tion [3] below, the resulting cascade decay to the ground 
state contributes the line emission. In a CIE plasma, 
on the other hand, origins of the line emissions are more 
dominated by the collisional excitation. Therefore, a kT z 
determination done by comparing with a CIE plasma 
code is not a proper method. Also, the previous works 
assumed that the continuum spectrum purely consists 
of bremsstrahlung emission, and determined kT e to be 
~ 1.0 keV. This value has been reduced owing to the 
discovery of the strong RRC emissions. 

It should be noted that the elemental abundances in 
Table[T]are determined as the parameters of the 0.6 keV 
VAPEC (CIE) component, and hence should be modified 
in the real case of the overionized plasma. The intensity 
of the Si-Hea line is given as Fn ea oc s(T e ) ■ fne{T z ) ■ Zsi, 
where s(T e ) and fn c {T z ) are, respectively, emissivity 
coefficient for electron temperature T e and fraction of 
He-like ion for ionization temperature T z . The abun- 
dance, therefore, can be modified by the fraction ratio 
of He-like ions, /h c (0.6 keV)//He(l-0 keV). According 
to the ionization population calculations by Mazzotta 
et al. (1998), the ion fractions of He-like, H-like, and 
fully-ionized Si (/h c , fa, fo) are estimated to be (0.83, 
0.15, 0.01) for kT z = 0.6 keV and (0.37, 0.43, 0.19) for 
kT z = 1.0 keV, respectively. Then, the modified Z§i is 
0.82 x (0.83/0.37) ~ 2.2 solar. For S, (/ Hc , fn, fo) are 
(0.91, 0.02, 0) for kT z = 0.6 keV and (0.59, 0.32, 0.06) 
for kT z = 1.2 keV. Thus, we similarly modified Zg to be 
1.7 x (0.91/0.59) ~ 2.6 solar. 

Using the measured RRC flux, we can indepen- 
dently determine the volume emission measure, VEM = 
/ n e n p dV/(4:TrD 2 ), from a following equation: 



^rrc = ai(T e ) • n z /n p ■ f • VEM , 



(2) 



where a\{T e ) and nz are K-shell RRC rate coefficient for 
electron temperature T e and number density of element 
Z, respectively. According to Badnell (2006), we find 
total radiative recombination rate for fully-ionized Si, at 
kT e = 0.6 keV, is a tot ~ 2.3 x 10~ 12 cm 3 s" 1 . Since 
the rate of the recombination into a level of principal 
quantum number n is described as; 



a n oc 



1 /3fcT e 
2~J7 



1 



(3) 



(e.g., Nakayama & Masai 2001), we obtain cti/atot = 
0.65, for kT e = 0.6 keV and I z = 2.666 keV. Thus, the 
VEM is calculated to be ~ 9.9 x 10 12 (Z si /2.2 solar)" 1 
(Zo/0.19)- 1 cm" 5 . Similarly, a to t at 0.6 keV for fully- 
ionized S is derived to be ~ 3.2 x 10 -12 cm 3 s _1 . There- 
fore, the RRC flux of S corresponds to the VEM of 
~ 9.4 x 10 12 (Zg/2.6 solar)- 1 (/b/0.06)- 1 cm" 5 . These 



values are almost consistent with that of the VAPEC 
component (Table [1]). 

The rectangular region in Figure [T] corresponds to 
3.1 x 2.2 pc 2 at D = 1.5 kpc. Assuming the plasma 
depth of 3 pc, the emitting volume is estimated to be 
6.1 x 10 56 cm 3 . Therefore, the VEM of 6.6 x 10 12 cm" 5 
is converted to the uniform electron density of n e ~ 
1.7 cm~ 3 . 

As the mechanism to form the overionized plasma, 
Kawasaki et al. (2002) proposed that the SNR consists 
of a central hot (kT e ~ 1.0 keV) region surrounded by 
a cool (kT e ~ 0.2 keV) outer shell, and interpreted that 
the hot interior cooled down via thermal conduction to 
the cool exterior. Under the several reasonable bound- 
ary conditions, they estimated the cooling time from kT e 
= 1.5 keV to 1.0 keV is about (3-10) x 10 3 yr, roughly 
the same as the SNR age (~4000 yr: Troja et al. 2008). 
However, this scenario cannot work on our new results. 
According to equation (5) of Kawasaki et al. (2002), the 
conduction timcscale is estimated to be 



t cond £Si 3.9 x 10' 



I'- 
ll. 7 



kn 



0.6 keV 



-5/2 



yr, (4) 



if the similar boundary condition is assumed. Thus, cool- 
ing via conduction requires far longer time than the SNR 
age. Photo-ionization is also unlikely because no strong 
ionizing source is found. Furthermore, the temperature 
of ^0.6 keV is significantly higher than that of a typical 
photo-ionization plasma (<0.1 keV: e.g., Kawashima & 
Kitamoto 1996). 

Since the progenitor of IC 443 has been suggested to be 
a massive star with strong stellar wind activity (Braun 
& Strom 1986; Meaburn et al. 1990), we propose another 
possibility that the rapid and drastic cooling is due to a 
rarefaction process, as discussed by Itoh & Masai (1989). 
If a supernova explodes in a dense circumstellar medium 
made in the progenitor's super giant phase, the gas is 
shock-heated to high temperature and significantly ion- 
ized at the initial phase of the SNR evolution. Subse- 
quent outbreak of the blast wave to a low-density ISM 
caused drastic adiabatic expansion of the shocked gas 
and resultant rapid cooling of the electrons. The life- 
times of the fully-stripped ions are roughly estimated to 
be t = (atot^e)^ 1 ^ 8.1 x 10 3 K/1.7cm~ 3 )- 1 yr for Si 
and ~ 5.8 x 10 3 (n e /1.7cm — 3 ) — 1 yr for S, respectively. 
Note that these values are underestimated compared to 
the actual timescale for the plasma to reach CIE, because 
the contribution of collisional ionization processes is ig- 
nored. Nevertheless, the estimated lifetimes are longer 
than the age of IC 443. The overionized plasma can, 
therefore, still survive at present. 

Future observations with very high energy resolution 
like the Astro-H mission will give firm evidence for the 
RRC structure not only on Si and S but also on the 
other major elements. This will provide more quantita- 
tive study on the peculiar SNR IC 443. 
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